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We study the complex visibilities of the cosmic microwave background 
anisotropies that are observables in interferometric observations of the cosmic 
microwave background, using the multipole expansion methods commonly 
adopted in analyzing single-dish experiments. This allows us to recover the 
properties of the visibilities that is obscured in the flat-sky approximation. 

t~^ \ Discussions of the window function, multipole resolution, instrumental noise, 

^■j I pixelization, and polarization are given. 
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I. INTRODUCTION 

The detection of the large-angle temperature anisotropy of the Cosmic Microwave Back- 
ground (CMB) by the COBEDMR experiment |l|] provided important evidence of large-scale 
spacetime inhomogenities. Since then, many CMB measurements have reported detections 
or upper limits of the CMB anisotropy power spectrum over a wide range of scales p[. Re- 
cently, BOOMERANG || and MAXIMA [|] data have revealed the structure of the first 
Doppler peak which arises from acoustic oscillations of the baryon-photon plasma on the 
last scattering surface. Future ground-based and balloon-borne experiments, and upcoming 
NASA MAP satellite will unveil detailed features of the CMB anisotropy, thus allowing one 
to determine to a high precision a number of cosmological parameters . 

The CMB is linearly polarized when the anisotropic radiation is scattered with free 
electrons near the last scattering surface 0. The degree of polarization is about one-tenth 
of the temperature anisotropy at sub-degree scales j^. The CMB polarization contains a 
wealth of information about the early Universe as well. It provides a sensitive test of the 
reionization history as well as the presence of non-scalar metric perturbations, and improves 
the accuracy in determining the cosmological parameters |^. So far, the current upper limit 
on the CMB linear polarization is 16 fiK 0. A handful of new experiments, adopting low- 
noise receivers as well as long integration time per pixel, are underway or being planned [ITU 



CMB experiments are commonly single-dish chopping instruments, whose scanning strat- 
egy and data analysis procedure are well developed. In the last decade, interferometers were 
introduced to study the microwave sky. Recent advancement in low-noise, broadband, GHz 
amplifiers, in addition to mature synthesis imaging techniques, has made interferometry a 
particularly attractive technique for detecting CMB anisotropics. An interferometric array 
is intrinsically a high-resolution instrument well-suited for observing small-scale intensity 
fluctuations, while being flexible in coverage of a wide range of angular scales with the res- 
olution and sensitivity determined by the aperture of each element of the array and the 
baselines formed by the array elements. A desirable feature of the interferometer for CMB 
observation is that it directly measures the power spectrum of the sky. In addition, many 
systematic problems that are inherent in single-dish experiments, such as ground and near 
field atmospheric pickup, and spurious polarization signal, can be reduced or avoided in 
interferometry. A brief account of the interferometic CMB observations can be found in 
Ref. |jll|. Currently, several new interferometers including the Very Small Array (VSA) ||12|| , 



Degree Angular Scale Interferometer (DASI) |jT3|, Cosmic Background Imager (CBI) [[14 



Array for Microwave Background Anisotropy (AMiBA) |T5| are being built, with sensitivity 
about a few fiK per beam. As for the polarization capability, the AMiBA will have full 
polarizations, whereas the VSA, DASI, and CBI will not be polarization sensitive initially. 
The output of an interferometer is the visibility that is the Fourier transform of the 
intensity fluctuations on the sky. In this paper, we will study the CMB anisotropics in 
the visibility space. There are several papers dealing with CMB anisotropy data from 
interferometers [^^-|T^ , and making explicit contact of the visibility with the angular power 



spectrum in /-space that is frequently used in single-dish experiments [p^ , |20| -p^] . All of 
them have treated the CMB sky as fiat due to the fact that typically a small field on the 
sky is viewed by the interferometer. This turns the analysis into a familiar two-dimensional 
Fourier transform problem. Here rather than assuming the fiat-sky approximation, we will 



perform a full-sky analysis of the visibility. Although we then have to give up the relatively 
simple Fourier-transform formalism, the bonus is that the angular power spectrum can be 
directly transferred to the visibility space. Hence, the statistics of the CMB visibilities is 
straightforward induced by the Gaussian variables in /-space. 

The paper is organized as follows. In Sec. ^ a brief account of the building block of CMB 



anisotropies is given. We introduce the interferometric observation in Sec. |T|, and present 
the calculations of the CMB visibilities in Sec. |V|. The resolution in /-space is discussed 
in Sec. 0, the instrumental noise and pixelization is introduced in Sec. ^I|, and the CMB 



polarization is treated in Sec. |V1J| . Estimation of signal to noise ratios for interferometric 
CMB measurements are made in Sec. IVllll. Sec. (TXI is our conclusions and discussions. 



II. CMB ANISOTROPIES 

Before we study the interferometric observation of CMB anisotropies, let us review the 
basic ingredients that are based upon the multipole expansion of the Gaussian random fields 
and commonly used in analyzing single- dish experiments 

Polarized emission is conventionally described in terms of the four Stokes parameters 
(/, Q, U, V), where I is the intensity, Q and U represent the linear polarization, and V de- 
scribes the circular polarization. Since circular polarization cannot be generated by Thomson 
scattering alone, the parameter V decouples from the other components and will not be con- 
sidered. Let us define T be the temperature fluctuation about the mean, then the CMB 
anisotropies are completely described by (T,Q,U), where each parameter is a function of 
the pointing direction e on the celestial sphere. 

Considering the CMB as Gaussian random fields, we can expand the Stokes parameters 

11,11] 



as 



Tie) = ^aT,iniYimie), 

Im 

{Q -iU){e) = Xl«2,«m2^«m(e), 
Im 

{Q + tU){e) = J2a-2M-2Yim{e), (1) 



Im 



where aT,im and a±2,im are Gaussian random variables, and ±2Yim are spin-2 spherical har- 
monics. Details of the spin-weighted spherical harmonics and their properties can be found 
in Appendix A. 

Isotropy in the mean guarantees the following ensemble averages: 



'^T,l'm''^T,lm/ — CTlSl'lSm'm, 

2,l'm''^2,lm) ~ {^ El + C Bl)5in^m'm, 

yh.,l'm'(^-'i,lm) — {CeI — CBl)5lil5m'm, 

'^T,l'm''^2,lm) — ~CciSl'l5m'm, (2) 



where Cti, Ceu Cbi, and Cci are respectively the anisotropy, E-polarization, B-polarization, 
and T-E cross correlation angular power spectra. 



Consider two pointings ei and 62 on the celestial sphere. Using the generalized addition 
theorem (|73D and Eq. @, we obtain the correlation functions |25], 



2/ + 1 
(T*(ei)T(e2)) = E ——CnPiicosP), 



(/-2)! 



CCTP,'(cos/3)e 



-2Ja2 



([Q(ei) + ^f/(ei)]*[g(e2) + zf/(e2)]) = E 
([g(ei) - 2f/(ei)]*[Q(e2) + tUie^)]) = E 



/2Z + 1 
An 

I2I + I 
An 



{CEi + CBi)2Yi-2{P,0)e 



-2i(a2-ai) 



{CEi-CBi)2Yi2{P,0)e 



-2i{a2+ai) 



(3) 



where (3, ai, and 02 are the angles defined in Appendix A. Therefore, the statistics of the 
CMB anisotropy and polarization is fully described by the four independent power spectra 
or their corresponding correlation functions. The details about the evaluation of the power 



spectra can be found in Refs. p3| , pl |. 

In realistic CMB observations, due to the finite beam size of the antenna and the beam 
switching mechanism, a measurement is actually a convolution of the antenna response and 
the Stokes parameters. This can be accounted by a mapping of the spherical harmonics in 
Eq. ®, 

sYUe)^i-irsWisYiUe), (4) 

where sVF/m is the window function. For a simple single-dish experiment with Gaussian 



beamwidth 0"^, it was found that ^5\ 



,Wira ^ exp [- (/(/ + 1) - s^) al 



(5) 



III. INTERFEROMETRIC OBSERVATION 



In contrast to single-dish experiments that measure or differentiate the signals in indi- 
vidual dishes, an interferometer measures correlation of the signals from different pairs of 
the array elements. Let us consider a two-element instrument and a monochromatic elec- 
tromagnetic source. The output, called the complex visbility, of the interferometer is the 
time-averaged correlation of the electric field E measured by two antennae pointing in the 
same direction to the sky but at two different locations ||26||: 



V{u) = {EiE;) = / deA{e,eo)Iie)e 



J2-Kiu-e 



(6) 



where u is the separation vector (baseline) of the two antennae measured in units of the 
observation wavelength, A denotes the primary beam with the phase tracking center pointing 
along the unit vector cq, and / is the intensity of the source. Note that u is generally a 
three-dimensional vector, and that V*{u) = V{—u) since A and / are real functions. 



In the following, we consider a primary antenna with Gaussian beamwidth at given by 

A(e,0) = e^, (7) 

where {6, (p) are the polar angles with respect to cq. Antenna theory |^ states that 




Ae//fiA = A^ with nA= f dnA{9,^), (8) 

where A is the observation wavelength, Qa is the field of view, and the effective area is the 
aperature efficiency times the physical area of the antenna: 

Aeff = VaAphy. (9) 

Let us assume a circular dish with diameter D, then we have 

(10) 

In typical interferometric measurements, we have A ^ D, dictating a small Qa- Thus, 
for a single pointing, it is very good to make the fiat-sky approximation by decomposing 

e = Co + X, with x ■ cq = 0, and |x| ^ 1, (11) 

meaning that x is a two-dimensional vector lying in the plane of the sky. Hence, the complex 
visibilty is reduced to the two-dimensional Fourier transform of the sky intensity multiplied 
by the primary beam: 

y(u) = y(u)e-'™-^« - /'dxA(x)/(x)e'™-", (12) 

where u is the two-dimensional projection vector of u in the x-plane. 

In a single-dish experiment, the resolution of the image is limited by the finite pri- 
mary beamwidth. In contrast, the interferometric imaging has resolution of the synthesized 
beamwidth which is determined by the sampling of the visibility plane (u-plane), where the 
minimum spacing is the physical size of the dish while the maximum spacing is limited by 
the size of the platform which resides the dishes. The sampling can be described by a sam- 
pling function S{u), which is zero where no data have been taken. One can then perform 
the inverse Fourier transform to form the so-called dirty image, 

J^(x) = /'rfur(u)S(u)e'2™.x^ ^^3^ 

Using the convolution theorem, I^ = AI * B, where -B(x) is the synthesized beam that is 
the Fourier transform of S'(u). So, in the fiat-sky approximation the calculation is simplified 
into a two-dimensional Fourier transform problem. However, the CMB is a large source. 
Although it is valid to treat the sky as fiat due to the small field of view of the primary 
beam in a single pointing, a larger scale CMB image requires multiple pointings on the sky. 
In addition, after replacing locally the sphere by a plane, the global property of the CMB 
field is obscured although there is a direct link between the angular power spectrum and 
the power spectrum in the plane. Therefore, in the following we will begin with the general 
form for the complex visibility in Eq. (^. 



IV. COMPLEX VISIBILITY OF CMB ANISOTROPY 

The CMB brightness fluctuation is related to the temperature fluctuation by 

f) R f)R 

AB, = -^AT, with ^^99.27/(a;)Jysr-VK-\ (14) 

where B^, is the Planck function of the photon frequency u, and 

f(-) = J^^.^ where - - 1-76 (y^) ■ (15) 

We write AT = T and then do the anisotropy angular power spectrum expansion (|l|). 

A. Infinite Resolution Limit 

Let us concentrate on the power spectrum and neglect the effect of the primary beam 
by taking for now the inflnite resolution limit that A = 1, then from Eq. (|]) the complex 
visibilty of the CMB anisotropy is written as 

V{u) = —^Vt{u), where Vt{u) = / deT(e)e2™-^ (16) 

Also expanding the phase factor in terms of spherical waves, 

^.iu-e ^ 47rY^^^M27ru)Y;^{e)YUu), (17) 

Im 

where ji is the spherical Bessel junction and u = \u\, we obtain 

Vt{u) = 47r ^ aT,im^'ji{27ru)Yi^{u). (18) 

Im 



For a flxed length of the baseline, Eq. (piq ) is analogous to a measurement of the temperature 
anisotropy T(e) in real space with a window function. So, the statistics of Vr on a sphere 
of radius u in the visibility space is completely described by the aT,im and the positively 
deflnite function 

Ii{u) = 167r^jf{2nu). (19) 

We will call this function the interference function. Estimation of the C^j's can be made 
from the visibility data in the same way as one does for an anisotropy sky map, with the 
sample variances determined by the coverage of the visibility sphere. 
Using Eq. (Q), the two-point correlation function is given by 

{V^{ui)Vt{u2)) = ^ E(2^ + ^)CTih{u)Pi{^o^B). (20) 

47r ; 

where P; is the Legendre polynomial and cos 6* = u\ ■ ii^- Hence, the rms temperature 
anisotropy in a given visibility is 
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vims = {V*{u)Vt{u)) = ^ Y.^21 + \)CtiIi{u). (21) 



On the other hand, based on the flat-sky approximation, the authors in Ref. [^ have 



obtained a similar autocorrelation function but with a different interference function, 

ll\u) = -J^i+iiAnu), (22) 

u 

where J2/+1 is the Bessel function. A comparison between these two interference functions 
is made in Fig. |l] with 2?™ = 200. Both interference functions have a high peak at / ~ 200 
and drop rapidly as / > 200. While Ii is positively definite, /; is oscillatory about the zero 
and has a lower peak. As such the latter would underestimate the rms fluctuation in a given 
visibility. For instance, assuming a flat anisotropy power spectrum, we find that the flat-sky 



rms visibihty is about one-fourth of that in Eq. (pT]) . 



B. Finite Primary Beam 

In practice, the instrument is limited by a finite primary beamwidth given by the aperture 
function. Inserting the aperture function with tracking direction e in Eq. (0), we have 

Vt{u, e) = I de'A{e', e)T(e')e2™-^'. (23) 

Doing the expansions (|l|) and ([T7|) , assuming a Gaussian beam (|^), and using Eq. (|74|), 
Eq. (p3D becomes 



/ 47r / 477 
VT{u,e) =8n'^Y.Yl (^T,hm^i^^ih{'^'^u)Yi,^rn2{u)d ^^ +l V 2/ +1 ^ 

li'mi l2m2 ' -'- ' 2 1^ 

E -n.Y*m2ie) -mYir^Ae)wkhj2.m), (24) 

where the square root of the window function is 

_ 132 
wHhJ2,m) = r d(3 13 e~^^Y*^iP,0)Yi,U/3,0), (25) 





and sYim is the spin-s spherical harmonics (see Appendix A). Eq. ( p^ ) is the general re- 
sult for the complex visibility of the CMB observed by an interferometer with two identical 
Gaussian primary beams separated by a baseline u. However, it is difficult to get any useful 
information from its present form. In the following, we will pursue two configurations for 
which Eq. ( P^ can be boiled down to useful forms. The first one is the "close-packed" con- 
figuration with the dishes almost in touching, i.e., uX ~ D. This configuration is commonly 
adopted in CMB observations in order to maximize the size of the dish to obtain optimal 
sensitivity for signal detection. The second is the "widespread" configuration in which the 
length of the baseline is much longer than the diameter of the dish, i.e., uX 3> D. This 
configuration is needed when one wants to resolve the fine structure of an image. 



1. Close-packed 

Let us start by analyzing the interference function, which is composed of the spherical 
Bessel function ji{27iu). Typically, the baseline m :» 1. For large argument, ji has a sharp 
peak at / ~ 2ttu. The asymptotic expansion in the peaked region is given by p9 



jK2™)~ Vhl(/ + i)-Ui(-2^y), where 2™ = (/ + ^) + y(/ + ^)^ (26) 

and Ai{z) is the Airy function. Thus, for a fixed u, the location and width of the peak are 
approximately given by 

/pfc + 0.84 = 2™, ^^^ ^^/~4- (27) 

In addition, the height of the peak is given by 

Jij27ru) ~ V^2-^/;,Ui(-l); Ai(-l) ~ 0.54. (28) 

In Fig. |1|, we have used the asymptotic expansion (^) to reproduce the interference function 
in Eq. ([T9|) . It shows that the asymptote is a fairly good approximation. 

In the close-packed configuration, we denotes the minimum spacing Imin = 27rM, where 
from Eq. (0) 

D 2 rr 

/^i„~27r- = -W->l. (29) 

A ab'i r]a 



la 



Because the dominant contribution in the /2-summation in Eq. (plf) comes from I2 ^ Imin, 
we can approximate the associate Legendre polynomial in the integrand ( p5[ ) by pO[ , 

P,7™(cos/3)~/2""'J^(/2/3). (30) 

Hence, the integral ( ^5]) can be evaluated analytically for certain limiting values of li. The 
integration is detailed in Appendix B. Essentially, the window function can be approximated 

by 



WHh,l2,m) 



2/i + l(/i + m)! 
4:71 (/i-m)! 



2l2 + l{l2 + m)\ 
47r (/2-m)! 



=e 2 



{hkr+-2 v/27r 



|(/i-'2)V?_ (31) 



Note that the window function has a Gaussian width, Alw ~ ^min, and that it decreases 
with m rapidly. The latter enables us to retain only m = in the m-summation in Eq. (p4[), 
which then becomes 

VTi%e) ~ V2^a, J2 J2^TMm,Yl,^Ae)^'H2l2 + l)ji,i2nu)Pi,ie ■ ^)-i=e-^('^-'^)''^^ (32) 



In the close-packed configuration, we usually have the condition that Alw ^ ^h- As such, 
as compared to the window function, the interference function can be treated as a delta 

1 /3 

function at I2 = Ipk, where Ipk + 0.8/^^ = Imin- For Imin ^ 1, we have Ipk ~ Imin- Hence, we 
approximate 



%l2 



-.e 2 



^ai-'2)V2 



--(It -J ■ )^a^ 



Ir 



in Eq. (JS^) , which then turns into the final form, 



Vrie) ~ e" 



where the window function is 



,u-e ySyfafe 



5I^r,imy«m(e)Vrj 



1 
I 1 



Im 



Wi = e"('"^™'")^'^6. 



(33) 



(34) 



(35) 



The result (0) shows that the baseline vector appears only in the irrelevant phase factor, 
and thus that the close-packed interferometric beam is similar to a single-beam antenna 
undergoing chopping and wobbling. In fact, this similarity has been pointed out in an 
early CMB interferometric measurement [0. So, when we attempt to construct the two- 
point correlation function similar to Eq. (|20| ) with ImiI ~ \u2\ ~ D/X, we obtain a pure 
phase function which does not contain any useful information. The reason is simply that 
the correlation over the domain in the w-space spanned by Ui and U2, whose size is still 
comparable to the size of the dish, is almost a constant. In light of this, in Eq. (|3^) we have 
replaced Vt{u, e) by Vr(e) to refiect that it is more legitimate to deal with the visibility in 
e-space than in -u-space. Therefore, for the close-packed configuration with a fixed baseline, 
we should sample Vrie) at different parts of the sky by repointing the entire telescope, and 
analyze the data in the same way as in the single-dish observation. 



2. Widespread 



In this case, the dominant contribution in the /2-summation in Eq. ( plD comes from the 
range, Imin -C ^2 ~ 2'7iu, where we can approximate the associate Legendre polynomial in 
the integrand (^ by |g. 



P,™(cos/5) 



^r(/2 + m-H) 
vr r(/2 + f ) 



sin 2 p cos 



\k + l)P + im-^-)^ 



Hence, the window function (|25|) can be approximately evaluated as 

"2/1 + 1 r(/i + m + i)r(/i - m + 1) 



WHh,l2,m) 




'^^-\{h-h?^l _ 



(36) 



(37) 



We have detailed the integration in Appendix B. Note that the window function (^) has a 
very narrow Gaussian width A/iy ~ Imini when compared to the baseline length ~ 27rM. As 
such, we can make li ~ /2- Then, using the asymptotic formula 



r(/±m + l) 

r(/ + i) 



^+2 



±m- 



for / > 1, 



(38) 



we can further approximate the window function as 

wHh,l2,m) ~ (27r)-iafee"5('i-'2)'-.', (39) 

which is independent of m in the asymptotic hmit. Furthermore, when the condition Al^ <^ 
All is satisfied, the window function can be treated as a delta function, 

wHh,k,m)c:^{2'K)-la,6i,i,. (40) 

Plugging Eqs. (|40| ) and ( [76D in Eq. (0), we finally obtain 

Vt{u) ~ VSnab ^ aTM^^Jii'^^u)Yi^{u), (41) 

Im 

which does not depend on the tracking direction. The structure of this equation is similar 
to Eq. ([T8|) , and thus the discussions following Eq. ([T8|) in Sec. |1VA| are equally applied to 
Eq. (ilD. 



Although the visibility (|41|) is reduced by the finite beam size, it enables us to make 
independent measurements of the visibility by pointing the telescope to different uncorrelated 
patches of the sky. Then, we can extract the Cti^s from the visibility data on the u-sphere for 
each pointing. So, if we have A^ uncorrelated pointings and a full u-sphere coverage for each 
pointing, the cosmic variance of the estimation of Cti^s will be reduced by a factor of y/N. 
It is interesting to note that the inevitable cosmic variance in single-dish experiments due 
to the existence of a single universe can be in principle circumvented in the interferometry. 
However, the widespread condition requires (27ru)^/^ ^ 1, meaning a very high resolution 
scale at which the primary CMB anisotropies have very low power. 

V. INCREASING RESOLUTION 

We have learned in the previous section that the resolution which we have in /-space 
for a single pointing of the close-packed interferometer is equal to the size of the primary 
beam. However, we can increase the /-resolution by combining several contiguous pointings 
of the telescope. This is analogous to the Fraunhofer diffraction in optics, in which narrower 
interference fringes are obtained by using many apertures. In Ref. |T^, they have demon- 



strated a case in which the resulting aperture in the w-plane can be made much narrower by 
considering N x N pointings on a regular grid. Here we can show the increase of /-resolution 
by combining various pointings of Vt{u, e) in Eq. (|2^) . The simplest case is to sum over all 
pointing directions, 

J dnVriu, e) = SttV^ ^ aT,ln.^' Ji{2nu)YiUu) , (42) 

Im 



where we have used the orthonormality condition (p8| ) to evaluate the integration. This 
shows that the /-resolution is now determined by the interference function rather than the 
window function. It means that the resolution in / is increased from Alw to A//. Since 
Eq. ( ^21) is a general result, it can be applied to the widespread configuration in Sec. P^ V B 2| as 
well, where the /-resolution is already given by A//. Thus we conclude that in interferometry 
the /-resolution can be increased by combining telescope pointings up to a limit set by the 
intrinsic interference function. 
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VI. INSTRUMENTAL NOISE 

In the single-dish CMB experiment, a pixehzed map of the CMB smoothed with a Gaus- 
sian beam is created. In each pixel, the signal has a contribution from the CMB and from 
the instrumental noise. A convenient way of describing the amount of instrumental noise is 
to specify the rms noise per pixel apix-, which depends on the detector sensitivity s and the 
time spent observing each pixel tpix'- (Jpix = sj ^/t^. The noise in each pixel is uncorrelated 
with that in any other pixel, and is uncorrelated with the CMB component. Let ^pix be 
the sohd angle subtended by a pixel. Usually, given a total observing time, tpix is directly 
proportional to Q^pix- Thus, we can define a quantity w~^ , the inverse statistical weights per 
unit solid angle, to measure the experimental sensitivity independent of pixel size ||31||, 



w ^ = ftpixali^. (43) 

In interferometry, the instrumental noise is usually specified by the image sensitivity per 



synthesized beam area |28 



AJ = ^ ^^^^^^^ , ^ / ^ (44) 

VsVa Aphy ^2NbNp V^uUnt ' 

where ks is the Boltzmann constant, Tsys and rjs are respectively the system temperature 
and efficiency, Nh and A''^ are respectively the numbers of baselines and polarizations, A^, 
is the bandwidth of observation frequency, tint is the integration time, and rja and Ap^y 
are defined in Sec. pT[ For example, the number of baselines formed by Na antennae is 
Nt = KiNa - l)/2. In Eq. 



With regard to CMB observations, because a given baseline u is sensitive only to a narrow 
range of / centered at about 2Tru, it is more suitable to use the sensitivity per baseline per 
polarization given by Eq. (|4^ ) with A*";, = 1 and Ap = 1, 

1 2kBTsvs 1 / X 

Si, = f 'y' __ . (46) 

risTja Aphy V2A^ 

Let us consider a simple two-element interferometer with baseline u. For the one in 
close-packed configuration, it is convenient to describe the instrumental noise on a pixelized 
sky map as in the single-dish experiment, with Qpix = ^a and apix = Sb/^ytpix, for which 
w"^ is specified on the celestial sphere. In the widespread configuration, we can adopt the 
same strategy but now on a pixelized u-sphere. The size of the pixel can be chosen as the 
resolution in u, i.e, flpix ~ (Au/w)^, and the noise in this pixel is also apix = S];,! ^ft^. Since 
/ ~ 27rM, we have '^pix ~ (A////)^. Hence, we can assign w~^ on the u-sphere. 

VII. POLARIZATION 

It is a routine for radio interferometers to measure the polarization of the radiation 



field ||3^. Polarization measurements are generally made using a pair of feeds on each 
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antenna. Usually, the feeds are sensitive to orthogonal circular or linear polarizations. For 
instance, if the dual-polarization feeds measure the right and left circular polarizations, 
then the output will be the four correlations: (RR*), (RL*), (LR*), and (LL*). They can 
be related to the four Stokes parameters {I,Q,U,V). Denoting their associated visibilty 
functions by (V/, Vq, Vu, Vy), we have 

{RR*) = Vj + Vv, 

(LL*) = Vj- W, 

{RL*) = Vq + tVu = V+, 

{LR*) = Vq-zVu = V_, (47) 

where we have neglected the parallactic angle of the feed with respect to the sky and the 
leakage from one polarization channel to the other polarization channel. Similiar to Eq. (P^, 
we have 

V±{u) = fde'A{e',e){Q ± if/)(e')e2™-^'. (48) 

It is difficult to analyze these equations in general. We thus proceed the analysis using 
the flat-sky approximation, i.e., m ■ e ~ 0. Let us first set up a rectangular coordinate 
{e^, iyjiz), and choose 6 = 6^, then it can be shown that (see Appendix C) 

4^6^™-^' ~ (27rM)^e-'''"e2™-^', (49) 

4^6^™-^' ~ (27rM)^e2^^"e2™-'^', (50) 

where 6'^ = m ■ e^^ is the polar angle of u in the w-plane parallel to the iy — iz plane, and the 
basis {iy, Cz) is used to define Q and U . As such, Eq. (|^) becomes 



V+{u) ~ (27rM)-V'^" / t/e'A(e',e)(Q + «f/)(e')4^e2™"', 

V_{u) ~ (27rM)-^e-2^^" f de'A{e' ,e){Q - zf/)(e')^4e'™"'- (51) 

Using Eqs. (|^) and (|72D , and following the steps to reach Eq. (p^ , we obtain 



V^{u) ~ Svr^e^^^^" ^ ^ %2,hmi^'^Jz.(2™) ^2Yi,„,,{u 

limi l2m2 



Att Att 



2/i + 1 V 2/2 + 1 
Y: -nXilm^ie) .mYi,rnAe)WHhM,m). (52) 



where 



W-^{h,h,m)= r d(3 (3 e~^^ ^2YiU(3,0)^2Yi,m{f3,0). (53) 

It has been shown that the window function for polarization measurements in single-dish 
experiments is well approximated by that for anisotropy as long as / ^ 2 [^ (also see 
Eq. d^)). This result can also be applied to Eq. ([53|) . So, we can approximate the window 
function (|53| ) by Eq. (p5|), and hence the subsequent analyses are the same as we have treated 
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the anisotropy in the previous sections. The only differences are the overall phase factor 
containing 9^-, the spherical harmonics of spin-2, and that u lies in the w-plane. 
For example, in the widespread configuration, we have 



V^{u) - V8^at,e^^''-Y.^TVm^'jii'^^u) ^2Yi 



lm(U 



(54) 



Im 



Analogous to Eq. (||), by using Eqs. (||) and (|73|) , we can construct four independent two- 
point correlation functions from Eqs. (^) and (Q). They are 

2 

{V^iui)VTiu2)) ^ ^E(2^ + 1)CtiIi{u)Pi{cos9), 



0", 






121 + 1 



Atx 



{CEl + CBl)Il{u)2Yl.2{e,^), 






crt 



2/ + 1 



(^;(.On(.^.))-^e-^E^^.(,^ 



(/-2)! 



CciIi{u)P^icost 



(55) 



where ImiI = 1^2 1 = u, and the separation angle 9 = 9u^ — 9u2- Hence, the rms total 
polarization in a given visibility is 



vl^s = (yi{u)v^{u)) 



(Tt 



iTC' 



Y,{21 + 1){Cei + Cbi)Ii 



iW- 



(56) 



VIII. ESTIMATION OF SIGNAL TO NOISE RATIOS 



In the previous sections, we have presented the basic results that can allow us to lay out 
the strategy in CMB interferometric observation, and to deal with the observed data. First 
of all, it is useful to make some simple estimates of the CMB signal to noise ratios for the 
upcoming CMB interferometers. 

For a close-packed interferometer, the CMB anisotropy signal in a given visibility is given 
by St = {dBjdT)VTrms/V2, where from Eq. (H, 



V, 



Trms 



{VmVT{e))-^ 



1 1 (Tfc 

1-^ 



V2ln 



E(2/ + l)CTie 



I, J _ \2 2 



(57) 



while the noise limit A^ is given by A/ (|4^) . Hence, we estimate the signal to noise ratio per 
single pointing as 



^ ^ i.6.v^„^^/W (^Y (W fl22^) (^) (^ 



N 



VlOGHz/ Vhr/ V 



u 



1 1 



.(/iK)2 



Y.i'^l + l)CTie 



_(l_j _ -)2 2 

yt bfmn j ^ fj 



(58) 



''rmn I 
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where f{x) is given by Eq. ([T5|). This formula is also applied for the CMB polarization 
except replacing Cti by Cei + Cbi- 

To estimate the S/N ratios, we simply neglect the B-polarization power spectrum, and 
approximate the anisotropy and E-polarization spectra for / = 300 — 1000 respectively by 

/(/ + 1)Cti ^ 27r(AT)2, AT ~ 50/iK; 

/(/ + 1)Cei ^ 2'k{AE)\ AE ~ 5/iK. (59) 

Let us consider a close-packed interferometer with 19 dishes in the hexagonal configuration. 
Then, it has 171 baselines in total, and the number of shortest baselines is 42. For D = 1.2m, 
z/ = 30GHz, Nh = 42, and Np = 2, the minimum spacing l^in = 2vrM ~ 27iD/X ^ 754. Hence 
the anisotropy and E-polarization S/N ratios are respectively given by 

^~3. /^ (^^y(t_^y(]^](^L_] 

N - ^''^^''" llOGHz j [hrj [ T,y, ) \bQ^^K) ' 

If we switch to a higher frequency z/ = 90GHz while fixing /^m — 754, then the dish size 
should be reduced to D = 0.4m. For Nt = 42 and Np = 2, we would expect the same S/N 
ratios as for u = 30GHz. However, corrected for the Rayleigh- Jeans limit, we find that the 
prefactor in the St/N is reduced to 2.5, whereas in the Se/N the prefactor is 0.25. 



IX. CONCLUSIONS AND DISCUSSIONS 

We have presented a full-sky analysis of the monochromatic CMB complex visibilities. 
First of all, an exact expression for the sky power spectrum is obtained in Eq. (^TJ). It 
has an advantage over the flat-sky approximation for the interference function being pos- 
itively definite, whereas the flat-sky interference function is rapidly oscillating about the 
zero. In the latter, care must be taken in summing over / for a flat spectrum due to signifi- 



cant cancellations |Tl|. Moreover, we have found that the fiat-sky approximation generally 
underestimates the power spectrum. 

A full-sky expression in Eq. ( P^ for the CMB complex visibility is our main result. It 
serves as the basis for analyzing the /-resolution in a given visibility, especially when a large 
sky scanning is needed in order to obtain a high resolution in /-space. We have shown in 
Eq. (^) that the full-sky scanning can increase the /-resolution from maximal A/ ~ D/X 
to A/ ~ /3. One should further check whether the resolution increases linearly with the sky 
coverage. 



We have worked out two limiting cases of the visibility equation (|2^) in which the statis- 
tics of the visibility becomes transparent. Firstly, we have shown that the close-packed 
interferometer is functioning like a single-dish switching experiment. Therefore, on the issue 
of obtaining hign /-resolution, it is important to study and compare the efficiency of the 
the usual method of synthesis imaging of the sky against that of the aforementioned sky 
scanning method. Secondly, we have suggested for the widespread configuration that one 
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should analyse the visibility data on the u-sphere. We have also pointed out that the in- 
terferometry can in principle avoid the cosmic variance by obtaining different u-spheres via 
multiple pointings of the telescope to uncorrelated patches of the CMB sky. It is interesting 
to study how to implement this concept in practical situations. 

In this paper, we have performed the calculations assuming a monochromatic electro- 
magnetic source, and allowing a nonvanishing geometrical delay r = Xu- eo, which measures 
the elapsed time for the wavefront reaching one antenna and then the other. However, when 
observing with a finite bandwidth A,^, one usually correlates signals at two separate points 
on the same wavefront in order to obtain full fringe. This can be done by including within 
the interferometric system a computer-controlled phase delay to compensate for r |^6[. As 
a consequence, the interferometer response is 

Viu, t) = / V{u, v)e^^'^''dv. (61) 

Ji/o-A^/2 

This is a Fourier transform with conjugate variables v and r, and can be inverted to extract 
the desired V{u, z/). 

Finally, we remark that CMB interferometric observations are radically different from 
traditional radio interferometry. We certainly need more studies on several important issues 
such as observing strategy, /-space resolution and mosaicing, optimal estimation of the power 
spectra, point source and other foreground subtraction, and ground pickup removal. 
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APPENDIX A: SPIN- WEIGHTED SPHERICAL HARMONICS 

The spin-weighted spherical harmonics are related to the representation matrices of the 
three-dimensional rotation group. If we define a rotation i?(a, /?, 7) as being composed of a 
rotation a around e^, followed by (3 around the new e' and finally 7 around e^', the rotation 
matrix of R will be given by |33] 



DL,^(a, /5, 7) = J^^ sYUP. a) e-^^\ (62) 



An explicit expression of the spin-s spherical harmonics is [| [p^,p^ 



^In Ref. [33|, the sign (—1)"* is absent. We have added the sign in order to match the conventional 
definition for Y/^- 
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' l-s" 



X 



E 



2/ + l(/ + m)! (/-m)! 
\ r + s — m 



21 



sm 



ly-'-' cot^'-^'-"^ 



!)■ 



where 



inax(0, ?7i — s) < r < min(/ — s,l + m). 



(63) 



(64) 



Note that the ordinary spherical harmonics F/^ = o^Zm- Using the expression (p5|), one can 
show the symmetry, 



/^^^(a,/3,7) = ^m-s(7,/5,a) 
They have the conjugation relation and parity relation, 

= (-l)'"+^_,l^ 



s^hni^^,' 



my 1 Y)i 



(65) 



(66) 



,Fz^(7r - ^, + 7r) = (-1)' -sYUe, 0). 
They satisfy the orthonormality condition and completeness relation. 



(67) 



(68) 



Y^ sY.lie'^cl)') .l^™(0,0) = 5(0' - (j))6{cose' - COS0). 

Im 



(69) 



Therefore, a quantity r/ of spin-weight s defined on the sphere can be expanded in spin-s 
basis, 



r]{9,(j)) =J2VlmsYlmid,(j)), 
Im 



(70) 



where the expansion coefficients rjim are scalars. 

The raising and lowering operators, 8 and d, acting on 77 of spin-weight s, are defined 
byH 



dr] = — (sin( 
dr] = — (sin( 



d .8 



[smt^j r], 



[smU) r]. 



(71) 



When they act on the spin-s spherical harmonics, we have |33 



dsYlm = -[(/ + S){1 - S + 1)]^ s^lYlm. 



(72) 
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Using these raising and lowering operations, one can obtain the generahzed recursion rela- 
tion 1 2^, which allows one to construct easily the high-/ spin- weighted harmonics from the 
low-/ harmonics. 

From the rotation group multiplication law, one can derive the generalized addition 



theorem [25 



Es.Y, 



ImK'^l^Vl) S2^lm\P2,V2 



Y,r 



/2/ + 1 



(-l)^-^^_,,rz..(/3,«2)e-^^^"\ 



(73) 



where /3 is the separation angle between the two pointings (^i, 0i) and (O2, (P2) on the celestial 
sphere. Connecting them by a geodesic, ai and 0:2 are the angles between the geodesic and 
the longitudes passing through (6'i, 0i) and {62, ^2) respectively. When si = S2 = 0, Eq. (|73D 
reduces to the familiar addition theorem for spherical harmonics. Similarly, we have 



s^/m(6'l, 01; 



47r 



21 + 



^E sYim'i/3,a2)e'''''-m'YU02,(l>2). 



(74) 



Furthermore, using the symmetry relation (|0^ and the addition theorem ([73|), we obtain 



121 + 1 



An 



-V 



51-S2 



s,YisA02-euO)e''''>''e-'''^\ (75) 



As long as the two pointings are identical, this becomes 

,x 2/ + 1 



7 . -m^i-sj 



— m-' 1—S2 



Att 



■'S1S2- 



(76) 



APPENDIX B: WINDOW FUNCTION 

The window function W2[li^l2,m) in Eq. (^^ is 

w'^ = f^d(3 (3 e"^y,:„(/?,o)r,,„(/5,o) 



2/1 + 1 (/i - m) 
Att (/i + m) 



2/2 + 1 (/2 - m) 
At: {I2 + m) 



where 



dppe ^P,7(cos/3)P,7(cos/3), 



(77) 



(78) 



and we restrict m > 0. This has already included the case with m < 0, because 
W-2ihj2,m) = Wh{hJ,,-m). 

Since the Gaussian function in the integral ( fTHD has a width of a?, ~ /^J^, we approximate 



the associate Legendre polynomial as [^ 
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Pr {cos (3) 



^r{i + m + r 



sin 2 p cos 



/T In ^ / 1 \ ^ 



for / > k 



Pf"" {cos (3) r^r'^J^nil (3) for /-/, 



Pr{cos(3) ~ (-1)' 



(/ + m)! ((3 



for / -C I 



■m\{l — m)\ \2 J 
Then, the integral can be integrated analytically for certain limiting values of /i and I2 

For h ~ Imin and k ~ Imin, 



(79) 



(/i + m 



(/i — m 
(/i + m 



{h — m 
{h +'m 



{h + ni 



{h -fn 
{k + m 



{h - m 



' dp P e ^P;7"^(cos/?)P;;"(cos/?) 

-K 



Ipl^ Jo 



ill -m)\ (/2-m)!/f/ 
where Im is the modified Bessel function, which has the limiting form 

1 



IrrXx) 



V2 



e^ for X > 1. 



''nx 



When lil2<jl S> 1, we have 



(/i + m)!(/2 + m)! 



CTb 



{h - m)\ {I2 - m)\ [l^l^y 



n^ 



=e 2 



\{h-l2?<yl 



Hence this gives the result in Eq. (pT|). 
For Ix > Ivain and h > Imin, 






m + 1 



COS 



\h + l)(3+{m-l)^ 



COS 



ir(/i+m + l)r(/2 + m + l^ 



TT r(/i + 



r(/2 



(/2 + -)/3 + (m--)- 
dp e ^"b cos [(/i - h 



r(/i + m + l)r(/2 + m + l) ab _i(,^_ 

r(/i + |) 



,2^2 



r(/2 + 1) 72^ 



h?^ 



(80) 



(81) 



(82) 



(83) 



Hence this gives the result in Eq. (^). 

Furthermore, we have found that when I2 ^ Imin, I is subdominant for /i < Imin- When 
h ~ Imin, I is subdominant for both /i :^ Imin and /i <^ /mm- 

APPENDIX C: FLAT-SKY APPROXIMATION 



We are going to evaluate 



2'^!'^ 



^27riuei-e2 



(84) 



under the condition that 61-62 — 0. This condition corresponds to the flat-sky approximation 
when u = uii is the basehne vector and 62 is the telescope pointing direction. The reader 
may refer to Refs. p3| , p^ , p!T| for different approaches. 
Using Eqs. ([I^), (^, and (^), we obtain 



^g2™ei.e2 =47r^i'j^(2 



TCU] 



2/ + l(/ + 2)! 
47r {1-2)1 



Yi-2{ei ■e2,ai) 



^ , , ,2/ + l (/ + 2)! ., 



(85) 



where we have substituted ai ^ 7r/2 under the flat-sky approximation. From the recursion 
relation, 



Pr-^Hx) + 2(m + 1] 



X 



VT 



x^ 



Pr^\x) + (/ - m)(/ + m + l)Pr{x) = 0, 



when X ~ and m = —2, we have 

Pl{x)c:^-il + 2){l-l)Pr\x). 

As such, Eq. ( pS] ) can be approximated as 

^e"-^-'^-'' ^ 4TrY,i'ji{2nu)l{l + l)Y:^{ei)YiUi2) ■ 



(86) 



(87) 



(88) 



Im 



Applying the operator (^ to Eq. (|88|) and doing the approximation (^) again, we find that 
^^g2™e,.e. ^ e-2^^" ^2'jz(2™)/2(/ + 1)2(2/ + l)Pz(ei ■ 62), (89) 



where 6^ = 7r/2 — 02- Since ji is a sharply peaked function at / ~ 2ttu for 2?™ ^ 1, we 
can approximate the factor P{1 + 1)^ in Eq. (^) by {27tu)^ and then take it out of the 
/-summation. Hence we obtain Eq. (^91) . We can follow the same steps to derive Eq. (pO|). 



19 



[1] 

[2] 
[3] 
[4] 
[5 

[6] 
[7] 
[8' 
[9 
[10 

[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 

[19 

[20^ 
[21 



REFERENCES 

G. F. Smoot et al, Astrophys. J. 369, LI (1992). 
L. A. Page, |astro-ph/ 99 ITTgUj 



A. E. Lange et al, |astro-ph/0005004 . 



S. Hanany et al, |astro-ph/0005125 



G. Jungman, M. Kamionkowski, A. Kosowsky, and D. N. Spergel, Phys. Rev. D 54, 

1332 (1996). 

M. J. Rees, Astrophys. J. 153, LI (1968). 

J. R. Bond and G. Efstathiou, Astrophys. J. 285, L45 (1984). 

M. Zaldarriaga, D. N. Spergel, and U. Seljak, Astrophys. J. 488, 1 (1997). 

C. B. Netterfield et al, Astrophys. J. 474, L69 (1995). 

S. T. Staggs, J. O. Gundersen, and S. E. Church, |astro-ph/990406^ . 

M. White, J. E. Carlstrom, M. Dragovan, and W. L. Holzapfel, Astrophys. J. 514, 12 

(1999). 

See [http://www.nirao.cani.ac.uk/telescopes/cat/vsa.htm] . 



See [http://astro.uchicago.edu/dasi/ 



See [littp:/ /astro. caltech.edu/^tjp/CBiy 



See [littp://www.asiaa.sinica.edu.tw/amiba/[ 

H. M. Martin and R. B. Partridge, Astrophys. J. 324, 794 (1988). 

P. T. Timbie and D. T. Wilkinson, Astrophys. J. 353, 140 (1990). 

R. Subrahmanyan, R. D. Ekers, M. Sinclair, and J. Silk, Mon. Not. R. Astron. Soc. 

263, 416 (1993). 

M. P. Hobson, A. N. Lasenby, and M. Jones, Mon. Not. R. Astron. Soc. 275, 863 (1995). 

M. P. Hobson and J. Magueijo, Mon. Not. R. Astron. Soc. 283, 1133 (1996). 

R. Subrahmanyan, M. J. Kesteven, R. D. Ekers, M. Sinclair, and J. Silk, Mon. Not. R. 

Astron. Soc. 298, 1189 (1998). 

[22] M. White, J. E. Carlstrom, M. Dragovan, and W. L. Holzapfel, [astro-ph/9912422 . 

[23] M. Zaldarriaga and U. Seljak, Phys. Rev. D 55, 1830 (1997). 

[24] M. Kamionkowski, A. Kosowsky, and A. Stebbins, Phys. Rev. D 55, 7368 (1997). 

[25] K.-W. Ng and G.-C. Liu, Int. J. Mod. Phys. D 8, 61 (1999). 

[26] A. R. Thompson, J. M. Moran, and G. W. Swenson, Inter ferometry and Synthesis in 

Radio Astronomy (John Wiley & Sons, New York, 1986). 
[27] J. D. Kraus, Radio Astronomy, 2nd ed., Cygnus-Quasar Books (Powell, Ohio, 1986). 
[28] J. M. Wrobel and R. C. Walker, Lecture 9, Synthesis Imaging in Radio Astronomy II, 

edited by G. B. Taylor, C. L. Carilla, and R. A. Perley, ASP Conference Series Vol. 180 

(Astronomical Society of the Pacific, California, 1999). 
[29] M. Abramowitz and L A. Stegun, Handbook of Mathematical Functions (1964). 
[30] L S. Gradshteyn and L M. Ryzhik, Table of Integrals, Series, and Products (Academic 

Press, 1980). 
[31] L. Knox, Phys. Rev. D 52, 4307 (1995). 
[32] W. D. Cotton, Lecture 6 of Ref. [g|]. 
[33] E. Newman and R. Penrose, J. Math. Phys. 7, 863 (1966); J. N. Goldberg et al., ibid. 

8, 2155 (1967). 
[34] R. Penrose and W. Rindler, Spinors and Space-time, Chapter 4 (Cambridge Univ. Press, 

1984). 

20 



FIGURE CAPTIONS 



Fig. |l| Solid and dotted curves are the interference functions Ii{u) with 2nu = 200 plotted 
using Eq. (|19|) and the asymptotic expansion p6| ) respectively. Dashed curve is the fiat-sky 
approximation in Eq. (B2). 
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